1 



Direct Detection of Vector Dark Matter 

Junji HiSANO^'^, Koji Ishiwata^ Natsumi Nagata-*^'^ and Masato Yamanaka^ 

^Department of Physics, Nagoya University, Nagoya 464-^602, Japan 
'^Institute for the Physics and Mathematics of the Universe, University of Tokyo, 

Kashiwa 277-8568, Japan 
_ ^California Institute of Technology, Pasadena, CA 91125, USA 

^ ■ '^Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 

^MISC, Kyoto Sangyo University, Kyoto 603-8555, Japan 

O 

^3 ■ In this paper, we complete formulae for the elastic scattering cross section of general 

vector dark matter with nucleons in the direct detection at the leading order of the strong 
coupling constant as, assuming that the dark matter is composed of vector particles and 
interacts with heavy fermions with color charge as well as standard-model quark. As an 
, application of our formulae, the direct detection of the first Kaluza-Klein photon in the 

^ l' minimal universal extra dimension model is discussed. It is found that the scattering cross 

I ' section is larger than those in the previous works by up to a factor of ten. 

(— I ' §1. Introduction 

\ Cosmological observations have established the existence of dark matter. Nowa- 

^> ' days, the energy density of dark matter is precisely determined using the Wilkinson 

. Microwave Anisotropy Probe (WMAP) satellite. The standard model of particle 

\ physics, however, can explain neither its existence nor its nature, which has been a 

lO ' mystery in particle physics and cosmology. The Weakly Interacting Massive Parti- 

. cles (WIMPs) in models beyond the standard model are a good candidate for dark 

matter. The relic abundance is naturally consistent with the cosmological observa- 
tion if WIMPs have TeV-scale mass and they are thermally produced in the early 
universe. This is the so-called thermal relic scenario. In this scenario, the dark 
matter is originally nonrelativistic and acts as cold dark matter in the era of the 
^ ■ structure formation of the universe. The leading candidate of WIMP dark matter 

is the lightest neutralino in the minimal supersymmetric standard model, which is 
a Majorana fermion. In addition, other models beyond the standard model at the 
TeV scale predict the existence of stable vector particles, such as the Kaluza-Klein 
(KK) photon^)' ^) in the universal extra dimensions (UED)^)'^^ and the T-odd heavy 
photon^^"^^ in the Littlest Higgs model with T-parity.^^' ^"^^ 

Now, the WIMP dark matter scenario has been tested in collider and direct 
detection experiments. In the Large Hadron Collider (LHC), (pair) production 
of the WIMP dark matter with TeV-scale mass is expected. On the other hand, 
XENONIOO,^^) which is the largest-volume detector ever, is in operation to detect 
it as scattering signals with nuclei on the earth. If the WIMP-scattering events are 
discovered and its properties are measured in the direct detection experiments, the 
measurements would be tested using the data that indicates dark matter production 
at the LHC. 



typeset using PTPTeX.cIs (Ver.0.9) 



2 



J. Hisano, K. Ishiwata, N. Nagata and M. Yamanaka 



With such recent progress on the experimental side, the theoretical prediction 
of dark matter signals must have better accuracy. The collider signatures of dark 
matter (such as in supcrsymmetric or UED models) have been intensively studied. 
On the other hand, the scattering cross section of dark matter with nuclei in the 
direct detection experiments has also been calculated more accurately. In recent 
works,^^)'-*^^) it was pointed out that a gluon-WIMP effective interaction is one of the 
leading contributors to the elastic scattering cross section, and it was also shown how 
to evaluate it precisely. (In those works, the results are shown for supersymmetric 
dark matter.) For the other models, in contrast, the gluon contribution and other 
leading terms arc not taken into account correctly in the calculation of the scattering 
cross section {^^' for the UED model and^^ for the Littlest Higgs model with 

T-parity). 

In this paper, we assume that WIMP dark matter is composed of vector particles, 
and evaluate its elastic scattering cross section with nucleons, which is relevant to 
the direct detection experiments. Here, wc consider the case in which vector dark 
matter interacts with standard-model quarks and exotic heavy fermions, which is a 
fundamental representation of SU{3)c- As an application of our formulae, we discuss 
the direct detection of the first KK photon in the minimal universal extra dimension 
(MUED) model. We found that the spin-independent scattering cross section with 
a proton is increased by up to a factor of ten, compared with those in the previous 
works. We also show that the spin-independent scattering cross section of the KK 
photon with a proton ranges from about 3 x 10~^^cm^ to 5 x 10~^^cm^ on the 
parameter region with the thermal relic abundance consistent with the cosmological 
observations. Future direct detection experiments with ton-scale detectors might 
cover the range. 

This paper is organized as follows. In the next section, we briefly review the for- 
mulation for the scattering cross section of general vector dark matter with nuclei. 
Next, in §3, we derive the effective Lagrangian describing the interaction of general 
vector dark matter with quarks and gluons in target nuclei. In §4, we review the 
MUED model and summarize the mass spectra of KK particles, and apply the formu- 
lae that are generally derived to the case of the first KK photon dark matter. Then, 
we numerically calculate the scattering cross section of the first KK photon dark 
matter with/without higher KK mode contributions, and discuss the significance of 
each effective operator. Section 5 is devoted to the conclusion. 

§2. Formalism of cross section for vector dark matter 

In this section, we formulate the elastic cross section of general vector dark 
matter with nuclei. The cross section is calculated in terms of effective coupling 
constants, which are given by coefficients of eflFective interactions of vector dark 
matter with light quarks {q = u, d, s) and gluons. To begin with, we write down all 
the effective interactions that are relevant to the scattering. In the derivation, we 
referred to reference.^®-' 

Since the scattering process is nonrelativistic, all the terms that depend on 
the velocities of dark matter particles and nuclei are subdominant in the velocity 
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expansion. Therefore, we neglect the operators suppressed by the velocities of the 
dark matter particles or nuclei in our study. Then, physical degrees of freedom of the 
vector field, which we denote as with the mass M, where we also assume a real 
field, are restricted to their spatial components since d^B^ = 0. As a consequence, 
in the expansion of a strong coupling constant and in the nonrelativistic limit of the 
scattering process, the leading interaction of the vector field with quarks and gluons 
is given as 

^eff^ ^ £f + (2-1) 

q=u,d,s 

with 

< = f^m.B'^B^qq + ffB^B^^qi^q + ^e^^p^B^^ EPq^l^'q + ■^B'^idHd'' B^Of,,, 

(2-2) 

£f = faBPBpG'^^^Gl, + j^BPid'^id'^BpOf,,, (2-3) 

where ruq is the mass of a quark and e^^'"^ is the totally antisymmetric tensor defined 
as e''"^^^ = +1. Here, the covariant derivative is defined as = + iQgA^Ta, with 
Qs, Ta, and A'^ being the SU{3)c coupling constant, generator, and the gluon field, 
respectively. Ofiu and Ofiu are the twist-2 operators (traceless parts of the energy- 
momentum tensor) for quarks and gluons, respectively, 

= \qi {d^^. + D,^^ - ^g^^p^ q , (2-4) 

Of,, ^ (G';fG% + i^^.GS^G'^-^) . (2-5) 

Here, G^^ is the field strength tensor of a gluon. In the derivation of the effective 
Lagrangian, we used the equations of motion for vector field, i.e., {d^ — M'^)B^ = 0. 
The third term in the right-hand side of Eq. (2-2) contributes to the spin-dependent 
interaction in the scattering of dark matter with nuclei, while the other terms in the 
effective Lagrangian contribute to the spin-independent one. 

Scattering amplitude is given by the matrix element of the effective interaction 
put between initial and final states. In this evaluation, we use the equation of motion 
for light quarks. The validity of the application is proven in the case where the 
operator is evaluated using on-shell hadron states. Then, denoting | A^) (N = p, n) 
as the on-shell nucleon state, we can evaluate the second term in Eq. (2-2) as 

{N\qipq\N) = {N\mqqq\N). (2-6) 

Now, let us first examine the spin-independent scattering. As in the scattering 
process, the momentum transfer is negligible; the matrix elements between initial 
and final nucleon states with the mass tun {N = p,n) are obtained as 

{N\mgqq\N)/mN = hq , 

^ - hq ^ JtG , 

u,d,s 
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{Nip)\OlJNip)) = -^{Pf^Pu - \mlg^,) {q{2) + q{2)) , 

{N{p)\0%\N{p)) = —{p,p. - \mlg^,) G{2) . (2-7) 

In the matrix elements of twist-2 operators, (?(2), q{2), and G{2) are the second 
moments of the parton distribution functions (PDFs) of a quark q{x), antiquark 
q{x) and gluon g{x), respectively, 

q{2) + ^(2) = f dxx [q{x) + q{x)] , 
1 

G(2) = [ dx X g{x) . (2-8) 
Jo 

Then, after the trivial calculation of matrix elements for vector dark matter states, 
we obtain a spin-independent effective coupling of vector dark matter with nucleons, 
/at, as 

/iv/miv= fifTo+ E ^('?(2) + ?(2))5, 

q=u,d,s q=u,d,s,c,b 

-^/tg/g + ^G(2)5g, (2-9) 

where fg = fjP' + fq^ and Og = g^/47r. Note that the effective scalar coupling of dark 
matter with gluons, fc, gives the leading contribution to the cross section even if it 
is suppressed by a one-loop factor compared with those of light quarks. On the other 
hand, the contributions with the twist-2 operators of gluons are of the next-leading 
order as go is suppressed by as. Thus, wc ignore them in this paper. 

The spin-dependent scattering can be dealt with in a similar manner; the spin- 
dependent effective coupling is given by 

ajv = dgAqN , (2-10) 

q=u,d,s 

with 

2s^AqN = {N\qj^j5q\N) , (2-11) 

where is the spin of a nucleon. 

For our calculation, we describe, in Table I, the input parameters for the hadronic 
matrix elements that we use in this article. The mass fractions and spin fractions (in 
a proton) of light quarks, frg and Aqp, are evaluated using the results in Refs.^^^'^'^^ 
and Ref.,^^^ respectively. The prescription for evaluating the mass fraction is de- 
scribed in these references and Rcf.^'^^ 

The second moments of PDFs of gluons, quarks, and antiquarks in a proton, 
on the other hand, have scale dependence. We use the values that are calculated 
at the scale n = mz {mz is Z boson mass). This is because all the terms (and 
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Table I. Parameters for quark and gluon matrix elements used in this paper, frq are evaluated 
using the results in Refs.^^'^^"' The spin fractions for a proton are obtained from Ref.^^^ The 
second moments of PDFs of gluons, quarks, and antiquarks in a proton are calculated at the 
scale ij, = mz using the CTEQ parton distribution.^^' 



For proton 


fru 


0.023 


frd 


0.032 


fTs 


0.020 


For neutron 


/tu 


0.017 


frd 


0.041 


hs 


0.020 



Spin fraction 


(for proton) 


Aup 


0.77 




-0.49 


Asp 


-0.15 



Second moment at /it = mz 

(for proton) 


G(2) 


0.48 






u(2) 


0.22 


u{2) 


0.034 


d{2) 


0.11 


d{2) 


0.036 


s{2) 


0.026 


s{2) 


0.026 


c(2) 


0.019 


c(2) 


0.019 


6(2) 


0.012 


5(2) 


0.012 



coefficients) in effective couplings should be evaluated at the scale of the mass of the 
particle that mediates the scattering.*^ It is obvious at tree level scattering. Even 
in the case where the loop diagram induces scattering, it is understandable from the 
fact that the loop integration is dominated by the particle mass.**) The numerical 
values of the second moments of PDFs are shown in Table I using the CTEQ parton 
distribution. The spin fractions and second moments for a neutron are given by 
the exchange of up and down quarks. 

Eventually, using /n and a^v, we express the elastic scattering cross section of 
vector dark matter with target nuclei as follows: 



a 



8 </ -|- 1 



(2-12) 



Here, tua is the mass of target nucleus A. rip and n„ are proton and neutron numbers 
in the target nucleus, respectively, while J is the total spin of the target nucleus and 
{Sn) = (^|5'jv|^) are the expectation values of the total spin of protons and neutrons 
in A. The first term in the brackets on the right-hand side comes from the spin- 
independent interactions, while the second one is generated by the spin-dependent 
one. 



*' Here, some readers may concern about it by considering (Q is momentum transfer) 

expansion in the operator-product expansion method, such as in the evaluation of the deep inelastic 
scattering cross section of hadrons. However, that is not what we are considering. Rather, we derive 
the effective Lagrangian in a similar manner to when one derives the weak-interaction Lagrangian 
at low energy from a fundamental theory. When one gives a theoretical prediction of physical 
observablcs from the weak-interaction effective Lagrangian, one should set a rcnormalization scale 
in the evaluation of the matrix element of the physical amplitude to the weak scale, i.e., the mass 
scale of the particle that mediates the scattering, in order to make the theoretical prediction accurate 
without calculating higher order quantum corrections. This is why PDFs at the weak scale are a 
good approximation. 

We thank Mark B. Wise for pointing this out. 
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B., 







Fig. 1. Tree level diagrams of color triplet fermion V'2 exchange to generate interaction of vector 
dark matter with quarks. 

§3. Effective Lagrangian 



In this section, we derive the effective Lagrangian for general vector dark matter, 
which was introduced in the previous section. As we described in the Introduction, 
we assume that vector dark matter has interactions with fermions with color charge, 
in which quarks are included. Denoting these fermions as ipi and -02 with masses mi 
and 1712, respectively, we start with a simply parametrized interaction Lagrangian of 
the form: 



C = ip2 (a^/^aViT^ + H2->Pi'l''l5)^iB^ + h.c. 



(3-1) 



Here, we do not consider the case in which vector dark matter couples to a pair of 
quarks, since otherwise it would become unstable. We take m2 > mi without loss 
of generality here and thereafter. We also assume, for simplicity, that tpi and "02 are 
color triplets. It is straightforward to extend to the other cases. 

First, let us examine the scattering with quarks. The vector dark matter is 
scattered by quarks at the tree level. The relevant interaction Lagrangian is given 
by taking tpi = q'm. Eq. (3-1), and the diagrams are shown in Fig. 1. After integrating 
out the heavy particle ■02! the coefficients in Eq. (2-2), i.e., scalar-, axial- vector-, and 
twist- 2-type couplings with quarks, are induced as follows: 



?)2 

%2q 



m2 



rrio 



M2 



( 2 



2 



M2)2 



ml 

t2( 



M2 ' 



9q 



(3-2) 
(3-3) 
(3-4) 



Here, we take the zero quark mass limit. As is obvious, the effective couplings 
obtained here are enhanced when the vector dark matter and the heavy colored 
fermion are degenerate in mass. 

Next, let us discuss the interaction with gluons. The vector dark matter does 
not couple with gluons at the tree level; however, it does at the loop level, as shown 



B„ 
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Fig. 2. One-loop contributions to scalar-type effective coupling with gluons. 



in Fig. 2. In those loop level processes, all the heavy particles '01 and 1^2 that couple 
with contribute, in addition to quarks. The contribution of each diagram, (A), 
(B), and (C), to the scalar-type effective coupling with gluons fc in Eq. (2-3) is 
written as 



(42V1 + ^%2^i)f+i^'^ mi,m2) + (a^2^^ - 6^2^J/P(M; mi, ma) 
where and /i^^ (I=A, B, and C) are given after the loop calculation, 



(3-5) 



/f\M;mi,m2) 



1 



2 



Qmlml[{ml - mjf - M'^^mf + 3ml)] 



mT 



12M4 \ml 



+ 



12Z\2M4 



ml + mi- M^)A^ + 2m2A{5m'^ + 20m(m^ - mf M^{9m^ + mf)} 



+ 12m2{M^(m2 + IQmim^ + 5mf) - {m^ - mifim^ + Sm^)} 



,"2 2 



2\2/ 2 



, (3-6) 



/i^)(M;mi,m2) 



m2 



6miZ\2 



Z\(2m2 +ml- 2M'^) + Gmfmliml -m{- M^) 
m\rn^{A -|- mf (m| — ml + M^)} 



+ 



Z\2 



L, (3-7) 



/f)(M;mi,m2) 



2M^ - 3(mf m|)M2 + (mf 



m: 



2\2 



6Z\M2 



+ 



-^In^ 



mt 



+ 



mn 



In 



mS 



12M4 \^m2y 12M4 \ml 
A{M'^ - mf - ml){ml + m^) + 4mfm|{(mf - m^)2 - 2]VP{mj + m^)} 

UAM^ 



L, (3-8) 
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and 



/f ^(M;mi,m2) = /f ^(M;m2,mi), 
/?)(M;mi,m2) = /i^^(M; m2, mi), 
fP{M;mi,m2) = 0. 



Here, we defined functions A and L as 



A{M] mi,m2) 
L{M;mi,m2) 



l\A\ 



2M 
In 



ml + 1712) + 

M2 + ^^ 



mi) , 



tan 



M2 



(ZA>0) 
(^<0) 



(3-9) 
(3-10) 
(3-11) 

(3-12) 
(3-13) 



The total scalar-type effective coupling with gluons is obtained by the summa- 
tion of /'^^•* for I = A, B and C. Here, one should be careful in dealing with the 
contribution of the diagrams in the case where ^1 is a quark. In such a case, there 
exist diagrams in which momentum with the scale of quark mass dominates the loop 
integral. We call them "long-distance" contributions. On the other hand, wc de- 
scribe the contributions from the diagrams in which the loop momentum with the 
scale of particle heavier than quarks dominates the integral as "short-distance" ones. 
When one derives the effective coupling in the three-flavour effective theory from the 
six-flavour full one, one should not include the contribution of light quarks in the 
long-distance diagrams, as discussed in Refs.^^^'^^^ To identify the long-distance 
contribution among diagrams, let us consider the case where V^i is a quark {i.e., 
mq <C m2). In this limit, the long-distance contribution in fc should approach to 
— -^fq, which is obtained by the calculation of triangle diagrams for emitting two 
gluons using the scalar-type effective coupling with quarks in Eq. (3-2);^^) 



/i^)(M;m„m2) 



/f^(M;ma,m2) 
/l^^(M;m„m2) 
/f)(M;m„m2) 



mn 



6(m2 -M2)2' 
m2 

~3mg(m2 -M2)' 



+ 



rrio 



Q{ml - M2)M2 ' 12M4 
mq{ml - 2M2) 
6m2(m2 — M2)2 ' 
1 1 



In 



2 1 



M 



2\ n 



6(m| - M2) 6M2 12M4 



In 



2 1- 



m: 



2 / J 



(3-14) 
(3-15) 

(3-16) 

(3-17) 

. (3-18) 



It is found that diagram (A) gives the long-distance contribution, whereas the others 
yield the short-distance ones. Therefore, fa should be given as 



Ja 



E 



V'2,q=all 



(3-19) 
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Here, Cg = 1 + llas{mq)/ Air is the additional long-distance QCD correction pointed 
out in Ref.^^) We evaluate it as Cc = 1.32 and Cb = 1.19 for as{mz) = 0.118, and 
also adopt ct = 1 for simplicity. On the other hand, if there are no quarks in the 
loop, all the diagrams indicate as short-distance contributions. In such a case, fc is 
simply obtained as 

fG= E (/cf^ + /g ^ + ^ (3-20) 

V'l,V'2 

§4. Cross section of Kaluza-Klein dark matter with nucleons 

In this section, we apply the above results to the evaluation for the direct detec- 
tion of the lightest Kaluza-Klein particle (LKP) in a model with UED.^^'^^ Among 
the variety of UED models, the MUED model, in which one extra dimension is 
compactified on an S^/Z2 orbifold and a dark matter candidate is included, is most 
extensively investigated.^)'^-' In the MUED model, the stabilization of LKP is en- 
sured by the KK parity conservation. It originates from the momentum conservation 
along the fifth-dimensional space. The KK parity —1 (-1-1) is assigned to particles 
with odd (even) KK number. This assignment and the KK parity conservation 
prevent the lightest KK-odd particle from decaying to the standard-model particles. 

The mass spectra of KK particles including radiative corrections were calculated 
in Ref.^) The work^) found that the LKP was the first KK photon in the broad range 
on the MUED-model parameter space, and would be a dark matter candidate. The 
relic abundance of the first KK photon dark matter was precisely calculated. ^^^^^^^ 
From those analyse, the compactification scale 1/R consistent with the cosmological 
observations is 600 GeV < 1/-R < 1400 GeV and the standard- model Higgs boson 
mass is rrih ^ 230 GeV. From the indirect low-energy experiments and corresponding 
analyses, which arc b ,S7,^^) B meson decay, B'^-B^ oscillation,^^) and elec- 
troweak precision measurements,'^'^) the lower bound on the compactification scale is 
obtained as 1/i? > 600 GeV. Thus, the favorable compactification scale in light of the 
relic abundance is free from experimental limits, and despite the simple framework, 
the MUED model is phenomenologically successful. 

In the following, we assume the first KK photon to be the LKP and evaluate 
its elastic scattering cross section with nucleons, which is relevant to the direct 
detection of the LKP dark matter. At the beginning of this section, we explain the 
parametrization of the MUED model, and describe the mass spectra of KK particles. 
Then, we give the effective couplings of the LKP with quarks and gluons. Finally, 
the numerical results are shown. 

4.1. Mass spectra of KK particles 

First, we summarize the mass spectra of KK particles relevant to the direct de- 
tection of LKP dark matter. We employ the minimal particle contents, i.e., gravitons 
and right-handed neutrinos are not included in this framework. Then, the first KK 
photon (7'-^)) or the first KK charged Higgs boson {H^^^")) becomes the LKP, de- 
pending on the compactification scale and standard-model Higgs boson mass. Since 



10 



J. Hisano, K. Ishiwata, N. Nagata and M. Yamanaka 



the H^^^^ dark matter scenario is excluded phenomenologically, we mainly study the 
parameter region where 7^^^ is the LKP unless otherwise mentioned. 

The mass of the first KK photon is obtained by diagonalizing the squared mass 
matrix for the first KK neutral gauge bosons, 



(4-1) 



This is written in a {B^^\ H^3(i)^ hasis, where B^^^i and VF^(^) are the first KK 
particles of U{1)y and SU{2)l gauge bosons, respectively. Here, gi (52) represents 
the U{1)y {SU{2)l) gauge coupling constant, and f(~ 246 GeV) is the vacuum 
expectation value of the standard-model Higgs boson. Radiative corrections, (5m^(i) 
and (^m^(i) , were calculated in the work^) as 

Here, A is the ultraviolet cutoff scale of the UED model, which is ordinarily taken 
to be yl 0{1Q/R)F> {((n) is Riemann's zeta function.) Since the off-diagonal 
component is so small in the mass matrix, the KK photon 7(1) is approximately 
identical to the weak eigenstate B^^\ Therefore, we consider B^^^ as the 7*-^\ and 
take the mass of the LKP equal to the mass of B^^\ On the other hand, the mass 
of the first KK charged Higgs boson is given by 

"i^±(i) = {^/Rf +m^w + S^HMi) ' (4-4) 

Sml,,,, = (^^gj + ^gl - A,) HA^R^). (4-5) 

Here, is the self-coupling constant of the standard-model Higgs boson, which is 
defined by Xh = mf^/v'^. When the standard-model Higgs boson is sufficiently heavy, 
the radiative correction Sm?jj^,^ gains a negative contribution so that the first KK 
charged Higgs boson is the LKP. 

Next, wc describe the mass spectra of KK quarks. They are important since 
most of the effective couplings of LKP dark matter with nucleons are very sensitive 
to the mass-square difference between the KK quarks and dark matter, which we will 
see later. The mass matrix of the n-th KK quarks is written in the weak eigenstate 
basis. The up-type quark mass matrix, for instance, is given by 

where u'^") and describe the n-th KK modes associated with the SU (2) ^ singlet 
up-type quark and the up-type quark in the SU{2)l doublet, respectively, and m„ 
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is the mass of the (zero- mode, i.e., standard- model) up-type quark. The down- 
type quark mass is written similarly. The parameters and are given as 
= n/R + SrUj^in) and m^^2 = n/R + 5m,jj(„) with the radiative corrections: 



n 



ln(yl2i?2 

R 167r2 



On the other hand, the radiative corrections in the mass matrix of the n-th KK 
down-type quarks {d'^'^^ and D'('*) for SU{2)l singlet and in doublet, respectively) 
are 

c A 2 ^ 1 2^ ^ ln(.l2i?2) 



4^V^ 167r2 ' 

5mjj(n) = 5mjj(n.). (4-8) 

The equality in the latter comes from the fact that L>'(") and make a SU{2)l 

doublet. In addition to these radiative corrections from gauge couplings, the third- 
generation KK quarks also receive the corrections from top Yukawa coupling yt. 
They are given by 

/ 3 2\nln(yl2i22) 

( ^ 2\n HA^R^) 

for the SU{2) singlet and doublet, respectively. The mass eigenstates of the n-th 
KK quark modes are expressed in linear combinations of their weak eigenstates. For 
example, in the up-type quark sector, the relationship between the mass eigenstates 
(denoted as n^") and ?/(")) and the weak eigenstates is expressed as follows: 



u^^^\ _ /— 75C0SQ;(") sina^^A / u*^") 



,?7'(")y V 75sinaH cos V^^^^^ 
Then, one obtains the mass eigenvalues from Eqs. (4-10) and (4-6) as 



(4-10) 



/ / (") I (") \ 2 (n) (n) 

with the mixing angle 

sin2aW= J"^" . (4-12) 
w-u + rnij 

The off-diagonal component in Eq. (4-6) is small so that the mass eigenstates are 
approximately given by the weak ones and the mixing angle is tiny. However, the 
mixing should be kept in finite when we calculate the scalar-type coupling constant 
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3. Mass degeneracy between LKP dark matter and first KK quarks, (mq(i) — M)/M, for 
A = 5/R, 20/R, and 50 /R. 



fq in Eq. (2-2), since such a type of coupling flips the chirality of the quark and the 
mixing is expected to contribute to the coupling. 

Figure 3 shows the mass degeneracy between LKP dark matter and first KK 
quarks, (m^(i) - M)/M, for q = U,u, T, and t. We set A = 5/R, 20/R, and 50/ R 
for panels (a), (b), and (c), respectively. In this calculation, the gauge coupling 
constants are set to be the values at the electroweak scale. For a larger cutoff 
scale, the mass degeneracy is relaxed owing to the logarithmic factor in radiative 
corrections. 

4.2. Scattering with light quarks and gluons 

Now, we calculate the scattering amplitude of the LKP with light quarks. This 

process is given at the tree level. Although several authors have already evaluated the 
contribution,^)' they mixed the contribution of twist- 2-type operators with 
the scalar-type one and did not evaluate them correctly. The relevant Lagrangian is 
obtained by taking (-i/-!, '02) = and {q,q'^^^) in Eq. (3-1). Here, g^") and Q(") 

describes the mass eigenstate of the n-th KK quarks, which are the SU{2)l singlet 
and doublet, respectively, and the coupling coefficients are given as (see Appendix) 

= -y(cosQ;(^)yqL + sinaWyqR), 
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Fig. 4. Tree level diagrams for elastic scattering of -S'^' with light quarks: (a) Higgs boson exchange 
contribution, and (b) KK quark exchange contributions. 



boWa = -^(-cosa«yqL + sina«yqR), 



91 



\Wq = y(-sma(i)yqL + cosa(i)yqR) 



(4-13) 



In the above expression, Y^jl and Y^r are hypercharges of left-handed and right- 
handed quarks, respectively. (For instance, Y^l = g5^R — ~\-) The mixing angle 
a^^^ is defined in Eq. (4-10). Those interaction terms generate diagrams of the first 
KK quark exchange (Fig. 4(b)). In addition, there also exists a tree level scattering 
process by the standard-model Higgs boson exchange (Fig. 4(a)). In the Higgs boson 
exchange process, not only the standard-model Higgs boson but n-th excited KK 
Higgs bosons (n > 2) propagate. However, this effect is so small compared with that 
of the standard-model Higgs boson that we ignore the contribution. 

Using the results in the previous section, we derive the coefficients of the effective 
Lagrangian in Eq. (2-2) as 



/. = - 



+ 



Ami 



9l 



m 



i^(m2 



'"qui 



igfM 



L^Q(l) 



gq = -gfM' 



-M2 



+ 



m 



(1) 



qR 



K(i) 



M2)2_ 



m 



(1) 



M2 



+ 



^qR 



(mJ(,,-M2)2 (mJ„-M2)2j- 



(4-14) 
(4-15) 
(4-16) 



The first term in Eq. (4-14) corresponds to the standard-model Higgs exchange con- 
tribution, while the other terms come from the KK quark exchange processes. 

Next, let us discuss the scattering of the LKP with gluons. The contribution 
of this scattering process has not been taken into account or has not been properly 
calculated before. As we will see, however, it is not negligible. This scattering pro- 
cess is induced by one-loop level. The diagrams that contribute to the scattering 
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9 9 9 9 9 

(A) (B) (C) 



Fig. 5. One-loop diagrams for the effective interaction of ' with gluons via KK quarks. Diagram 
(A) gives a long-distance contribution, whereas (B) and (C) give short-distance ones. 



with gluons arc categorized into three types: (i) quark and the first KK quark con- 
tribution, (ii) higher KK quark contribution, and (iii) standard-model Higgs boson 
exchange contribution, such that 

fG = f^+fP + ff^. (4-17) 

Let us first closely look at contributions (i) and (ii). They are induced from 
the diagrams depicted in Fig. 5. The type (i) contribution is given by the interac- 
tion Lagrangian, which was introduced in the previous tree level calculation (ie., 

(V*!) V'2) = (?) 9'-^^) and {q, Q^^^)). The relevant Lagrangian for the type (ii) contribu- 
tion, on the other hand, is given by taking (V'i,V'2) = (g''"'\ g*-"""*"^-*), (Q*^"'^ Q*^""*"^^), 
(g("),Q("+i)), and (Q("), g("+i)) (n > 1) in Eq. (3-1), with the couplings, 

a^(„+i)g(„) = - ^ [>"qL sin a(") sin a("+^) + FqR cos a^") cos a^^+^^J , 
v 2 

OqC'+DqW = -^[i^qLCOsa(")cosa("+i) +yqRsina(") 



&Q(n+i)g(n) = — ^[YgLsina^"-^ cosa^"'^^^ - YqECOsa^"^ sin 
V 2 



h^(r.+,)Q(u) = - [yqL COS Sin -YqR Sin COS 

6g(„+i)^(„) = 6q{„+i)q(„) = aQ(„+i)g(„) = a^{n+i)q{n) = 0. (4-18) 

(For the derivation, see Appendix.) 

As we discussed in §3, there exists the long-distance contribution when the 
quark runs in the loop. This is the case for the type (i) contribution. Thus, in 
the calculation of the type (i) contribution, the diagrams in which light quarks run 
in the loop should not be included. Then, the contribution to the scalar effective 
coupling to the gluon is written as 



.(i) _as sr 



q=c,b,t 



+ i'Ji),)/^ ^(M;m„m^a)) + (aj,,^ - 62^,)^)/i^)(M;m„m 
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Fig. 6. One-loop diagram for the effective interaction of B^^^ with gluons through the exchange of 
standard-model Higgs boson. 



g=all I=B,C 



i^QWq + ^Q(l)q)/+ (^;"^9'"^Q(l)) + (^QW? ~ ^Q(l)g)/- (-^; ^9' "^QW ) 



(4-19) 



On the other hand, for the type (ii) contribution, aU the diagrams are short-distance 
contributions. Then, 



.(11) _ as 



X] XI ] « [f+ (^i "^gC") ' ) + /P (M; m^(n) , m^(n+l) )] 

q=allI=A,B,C 



(4-20) 

Finally, we discuss the type (iii) contribution. The standard-model Higgs boson 
is coupled with the standard-model quarks and also higher KK quarks, and then the 
triangle loop diagrams of heavy quarks and KK quarks contribute to the effective 
scalar coupling to gluons. The corresponding diagram is shown in Fig. 6. In the 
triangle loop, all the higher KK quarks run in the loop, in addition to heavy quarks. 
Then, we obtain the standard-model Higgs exchange contribution as 



.(HI) ^ 



(cc + Cfe + ct) + mtct V sin 2a(") (^— + -^—) 

(Cc + Cb + Ct)+Ct'^ 



2ml 



(4-21) 



In this calculation, we ignore quark masses except that of the top quark. 
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4.3. Results 

Now, we are at the point of showing the numerical results of the cross section. 
At the beginning, let us examine the case in which n-th KK modes with n > 2 
are neglected so that we see features of the cross section. Also, it is worthwhile 
when one applies the formulae for other vector dark matter scenarios. In Fig. 7, the 
spin-independent cross section of the LKP with a proton (asi) is depicted. Here, 
we set m^(i) = mg(i) = mist for simplicity, except for the masses of the first KK 
top quarks, for which wc take m^(i) = ''T^'^ii) = "i^ist + "^t ("^t quark mass), 

and show the contours of the cross section on the plane, the LKP mass vs the mass 
degeneracy between the first excited KK quarks and the LKP, (mist — M)/M. The 
standard-model Higgs boson mass is set as m^ = 120 GeV, 200 GeV, and 500 GeV. 
Prom the figure, it is found that the cross section is enhanced when the first KK 
quark and LKP masses are degenerate, as is expected. Such a behavior was already 
discussed in the previous works. However, with the complete calculation that has 
been carried out here, we have discovered that the spin-independent cross section is 
larger than the former theoretical predictions. In the case of rrih = 120 GeV and the 
mass degeneracy of 10%, for instance, the cross section ranges from 3.5 x 10~^^ cm^ 
to 1.3 X 10^"^^ cm^ when the LKP mass is from 200 GeV to 1 TeV. This result is 
larger than those in the previous works^^'^^^'^^^ up to about a factor of ten. In our 
calculation, we have found that the twist-2-type operator contribution, which had 
not been taken into account correctly in the previous works, dominates the effective 
coupling when M < O(TeV) and enhances the cross section. In the parameter region 
M > O(TeV), on the other hand, the standard-model Higgs boson contribution (from 
both tree- and one-loop-level diagrams) dominates the effective coupling. Then, 
the cross section is determined using the standard-model Higgs boson mass in the 
parameter region. 

To clarify the behavior of the cross section, we also plot each contribution in 
the effective coupling /^v in Fig. 8. Here, we set nih = 120 GeV and 10% of the 
mass degeneracy. In the figure, "scalar" , "twist-2" , and "gluon" correspond to the 
first, second, and third terms in Eq. (2-9) (except for the standard-model Higgs 
exchange contribution), respectively, and the standard- model Higgs exchange con- 
tribution (including tree and one- loop levels) is denoted as "Higgs" . We have found 
that all the contributions have the same (negative) sign so that they are construc- 
tive. As a consequence, the cross section is enhanced in a wide parameter region. In 
the figure, it is seen that the twist-2 contribution dominates the effective coupling 
when M < O(TeV). The gluon contribution is sublcading, but not negligible at 
all. On the other hand, the standard-model Higgs boson contribution dominates the 
effective coupling when M > O(TeV). 

For completeness, we show the spin-dependent scattering cross section. The spin- 
dependent interaction is generated by KK quark exchange at the tree level, and it is 
found that the cross section is enhanced when the first KK quark and LKP masses 
are degenerate as expected from Eq. (3-3). In Fig. 9, the spin-dependent scattering 
cross section with a proton is shown as a function of the LKP mass assuming the 
degeneracies are 10%, 20%, and 25%. We found that the spin-dependent scattering 
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Fig. 7. Spin-independent cross section with a proton on the plane, LKP mass vs. mass degeneracy 
(mist -M)/M. Here, we take = 120 GeV (upper left), 200 GeV (upper right), and 500 GeV 



(bottom) . Lines correspond to the contours asi 
and 10"''® cm^ from bottom to top. 
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cross section is consistent with Refs.^-*'"^^)'^^-' 

Next, we discuss a more realistic situation for the direct detection of the LKP 
dark matter in the MUED model, including higher KK mode contributions. For 
the cutoff scale A = n/R, the KK particles up to the n-th mode are included in 
the calculations. Furthermore, we employ the radiative corrections for the spectra 
of KK particles discussed in §4.1. To see the importance of the higher KK mode 
contribution, we show the cutoff dependence of the spin-independent cross section 
with a proton by taking A = 5/R, 20/R, and 50/R in Fig. 10. In the figure, we 
take m/j = 120 GeV. For a larger cutoff scale, although a larger number of KK 
modes could contribute to the cross section, the cross section tends to be smaller 
than that for the smaller cutoff scale. This is because the large cutoff scale makes 
the mass degeneracy relax owing to the radiative corrections (as shown in Fig. 3), 
and this reductive effect dominates over the increasing effect from higher KK mode 
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200 300 500 1000 2000 

DM mass (GeV) 

Fig. 8. Each contribution in effective coupling fN/niN given in Eq. (2-9). Here, we set rrih = 
120 GeV and (mist — M)/M — 0.1. "scalar" , "twist- 2" , and "gluon" correspond to the first, 
second, and third terms in Eq. (2-9) (except for standard-model Higgs exchange contribution), 
respectively, and standard-model Higgs exchange contribution, including tree and one-loop lev- 
els, is denoted as "Higgs" . 



contributions. 

Finally, we discuss the feasibility of the direct detection of KK photon dark 
matter for the parameter space in which its relic abundance can be reproduced 
in accordance with the WMAP result. We plot the relationship between the spin- 
independent scattering cross section with a proton and the relic abundance in Fig. 11. 
The charged Higgs boson is the LKP in the upper right region, and hence this region 
is excluded. The dark (light) gray region shows the la {2a) allowed region from 
the WMAP observational result, which is the result calculated in Ref.^^-' Each line 
corresponds to the spin-independent cross section agi = 10~^^ cm^, 3 x 10~^^ cm^, 
10^^^ cm^, 4 X 10~^^ cm^, and 10"^'' cm^. The spin-independent cross section 
ranges from about 3 x 10~^^ cm^ to 5 x 10"^^ cm^ on the allowed region of the relic 
abundance, and are just below the present experimental bounds. The ongoing and 
future experiments for the direct detection of dark matter would reach sensitivities 
in this range. 

§5. Conclusions 

In this paper, we assume that WIMP dark matter is composed of vector particles 
and have evaluated an elastic cross section with nucleons in the direct detection 
experiments. The vector dark matter is predicted in models beyond the standard 
model, such as the KK photon in the UED model and the T-odd heavy photon in the 
Littlest Higgs model with T-parity. On the other hand, however, the cross section 
had not yet been consistently evaluated, even at the leading order of ag. We have 
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Fig. 9. Spin-dependent cross section with a proton. Here, we set (mist — M)/M = 0.1, 0.2, and 
0.25 from top to bottom. 




200 300 500 1000 2000 

DM mass (GeV) 

Fig. 10. Spin-independent cross section with a proton for nih = 120GeV. Each hne corresponds to 
A = ^/R, 20/ R, and 50/7?, respectively. 



derived the general formulae for the cross section of general vector dark matter with 
nucleons. 

As an application of our formulae, we discussed the direct detection of the first 
KK photon in the MUED model. We found that the cross section is larger than those 
in the previous works by up to a factor of 10. We showed that the spin-independent 
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Fig. 11. Contour plot of relic abundance consistent with WMAP result and spin-independent cross 
section in the (l/R, mn) plane for A = 20/ R. The dark (light) gray region shows the la 
(2a) allowed region in light of relic abundance. Each line corresponds to a spin-independent 
scattering cross section with a proton crsi = 10"''^ cm^, 3 x 10~*^ cm^, 10~*® cm^, 4 x 10"'*^ cm^, 
and 10""^ cm^ 



cross section of the KK photon with a proton ranges from nearly 3 x 10~^^ cm^ to 
5 X 10~^^ cm^ on the parameter region with the thermal relic abundance consistent 
with the cosmological observations. The future direct detection experiments with 
ton-scale detectors might cover the range. 

Note Added: While preparing the manuscript, we became aware of a paper 
by G. Belanger, M. Kakizaki and A. Pukhov.^^^ In that paper, the elastic cross 
section of the KK photon with nucleons in the MUED model was calculated using 

microOMEGA.^^-' In the calculation, the tree level contribution was included in the 
scattering process whereas the one-loop contributions were incomplete. 
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Appendix A 

Relevant Feynman rules 

In this Appendix, we collect parts of the Lagrangian in the MUED model. 
Standard-model Higgs boson coupling with the first KK particle of U (l)y gauge 
boson 

jObbho = Igfvh'BmBm^ (A-l) 
Standard-model Higgs boson coupling with n-th KK quarks 



■'q(iT')q{n)f 



n=l 



+ cos2a('^)(g('^)75Q(") - g("^759^"0 



(A-2) 



Coupling of the first KK particle of the U (l)y gauge boson with the first KK 
quarks and quarks 

^qM^B = -9iqi''[cosa^^^Y^LPL + sina^'^Y^KPRW^ Bj}^ 
-giQ'^'h^ [cos a^^^l^PL + sina(i)yqRPR]gS« 
+5197^ [sin a(i)yqLPL + cosa(i)FqRPR]g(i)S« 
+51^(1)7'^ [sin a(i)yqLPL + cosa(^)yqRPR]gS« (A-3) 

Coupling of the first KK particle of the U{1)y gauge boson with the n-th and 
(n -I- l)-th KK quarks 

Cg(n)gin+i)B = -^S«Xl[(^qLsina(") sina("+^) +FqRCOsa(") cosa("+^)) 

n=l 

x(g("+l)7^(^(«)+g(")7A'g("+l)) 

+ (FqL cos a(") COS a("+i) + YqR sin sin 

X (q("+i)7'^Q(") + Q("'7^Q("+^)) 
+ (FqLsina(") cosa("+i) - YqRCOsa^") sin 

X (Q("+^)7^75g^"^ + ^^"^7^75^^"+'^) 
+ (FqL COS a(") sin - YqR sin a^") cos a^^+i) ) 

X (g("+i)7'^75QW + 0(")7'^75g('^+^))] 

(A-4) 
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